Failure due to fatigue in fiber bundles and solids 
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Abstract 

We consider first a homogeneous fiber bundle model where all the fibers have got the same 
stress threshold (a c ) beyond which all fail simultaneously in absence of noise. At finite 
noise, the bundle acquires a fatigue behavior due to the noise-induced failure probability 
at any stress a. We solve this dynamics of failure analytically and show that the average 
failure time r of the bundle decreases exponentially as a — > a c from below and r = for 
a > a c . We also determine the avalanche size distribution during such failure and find a 
power law decay. We compare this fatigue behavior with that obtained phenomenologically 
for the nucleation of Griffith cracks. Next we study numerically the fatigue behavior of 
random fiber bundles having simple distributions of individual fiber strengths, at stress a 
less than the bundle's strength cT c (beyond which it fails instantly). The average failure 
time t is again seen to decrease exponentially as a — > cT c from below and the avalanche size 
distribution shows similar power law decay. These results are also in broad agreement with 
experimental observations on fatigue in solids. We believe, these observations regarding the 
failure time are useful for quantum breakdown phenomena in disordered systems. 
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I. Introduction 

If one puts a load or stress (a) on a solid or applies a voltage across an electrical circuit, a 
strain in the solid or a current through the circuit develops which grows linearly (Hooke's 
law or Ohm's law) with the stress or voltage. If the external load on the system increases 
beyond its threshold limit (<r c ), the system fails: stress a drops to zero due to fracture of the 
solid. The same occurs when the voltage on the network exceeds its limit and the current 
drops to zero due to the fuse of the circuit. Similar failures occur in dielectric materials when 
the electric field across the sample exceeds beyond its limit, and dielectric breakdown sets 
in. These failures usually nucleate around the defects in the solid and the failure behavior 
and its statistics therefore crucially depends on the disorder or impurity distribution within 
the sample. These (quasi-static) failure properties of disordered solids have been studied 
extensively in recent years [T]. 

The dynamics of these failures in such systems are quite intriguing and is being studied 
very intensively these days. The critical dynamics of failure and its universality class in the 
democratic (global load sharing) fiber bundle model [2] has been established very recently pj]. 
These dynamics of failure are intrinsic and induced by the successive stress redistributions 
due to the failure of weaker fibers. However, an important kind of dynamical failure due to 
fatigue |H occurs in such disordered systems when the fibers have an effective probability to 
fail under any stress |5J, or as the micro-cracks within the solid grow at the crack-tips with 
time due to chemical diffusion in the atmosphere The system then fails under a stress 
less than its normal strength (cr c ) and the time of failure (r) depends on the load applied on 
the sample: r 7^ for cr < a c and r ~ for a > a c . 

Here, we study first a phenomenological theory of crack nucleation, following Griffith 
IHE], at finite temperature (T) and estimate the average failure time r at any stress a less 
than a c . We then develop a simple model of fatigue-failure in a democratic fiber bundle 
model containing identical fibers of strength a c (homogeneous bundle), where the fibers 
have a finite noise-induced failure probability. We have derived analytically the failure time 
for the bundle as a function of the applied stress (a) and the noise (T). This result for the 
model is compared with that obtained for the phenomenological theory of crack nucleation at 
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finite temperature. It is also in broad agreement with some recent experimental observations 
on fatigue in disordered solids |UCj- Next, we derive the avalanche size distribution in this 
fixed strength model analytically and find robust power law decay. The above analytic results 
have been confirmed through the numerical studies on the same model. Finally we consider 
random fiber bundles with simple, yet nontrivial, distributions of the fiber strengths. Our 
numerical results show that for all these fiber bundles, the average time to failure r decreases 
exponentially as the stress level a approaches bundle's strength cT c from below and the 
avalanche size distributions show similar power law decay. We also discuss the plausibility of 
this (noise-induced) failure in other similar situations. In particular, we consider the validity 
of our model in quantum breakdown phenomena [8j: for example, in dielectric breakdown 
where the microscopic failure of the dielectric grains acquire a finite probability at any electric 
field due to quantum tunneling. The failure time and its variation with the strength of the 
external field in such a quantum failure can give us an estimate of the tunneling frequencies 
involved. 

II. Time for fracture in the Griffith nucleation model 

Griffith in 1920, equating the released elastic energy of a growing crack inside a solid with the 
energy of the newly created crack surfaces, came to a quantitative estimate of the fracture 
strength of a solid containing an already existing fixed geometry micro-crack. Assuming the 
linear elasticity behavior up to the breaking point of a brittle solid, the released elastic energy 
becomes E e i = (a 2 /Y)Iq for a three dimensional elastic solid under stress a, modulus of 
elasticity Y, containing a micro-crack of length l . The corresponding surface energy E s = §l\ 
where denotes the (crack) surface energy density. Using the concept of energy balance, 
Griffith equated the differential increment in the elastic energy dE e i with the corresponding 
surface energy increment dE s as the crack propagates a further length dl and got 



for equilibrium extension of the crack. Here a c is the amount of stress for and above which 
the micro-crack propagates in no time (or in a small time dependent on the sound velocity) 
and causes a macroscopic failure of the sample. 
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This quasi-static picture can be extended to fatigue behavior of crack propagation for 
a < cr c . At any stress a less than a c , the cracks can still nucleate [B] for a further extension 
at any finite temperature T with a probability ~ exp[— E/k B T] and consequently the sample 
fails within a failure time r given by 

r- 1 ~expi-E(l )/k B T], (2) 

where 

2 

E(k) = 4% - yll (3) 
is the crack (of length Iq) nucleation energy. Here k& is the Boltzman factor. One can 
therefore express (2) as 

2 

r~exp[A(l-^)], (4) 

where (the dimensionless parameter) A = '(YksT) and a c is given by (1). This immedi- 
ately suggests that the failure time r grows exponentially for a < a c and approaches infinity 
if the stress a is much less than a c when the temperature T is small, whereas r becomes 
vanishingly small as the stress a exceeds a c . 

III. Fatigue in a homogeneous fiber bundle 

Fatigue in fiber bundle model was first studied by Coleman in 1958 [2]. Thermally activated 
failures of fiber have recently been considered and approximate fatigue behavior has been 
studied We consider here a very simple fiber bundle model with noise-induced activated 
failure, for which the dynamics can be analytically solved. 

Let us consider a homogeneous bundle of iV fibers under load L{= A^cr), each having 
identical failure strength a c . Without any noise (T = 0), the model is trivial: the bundle 
does not fail (failure time r is infinity) for stress a < <j c and it fails immediately (r = 0) for 
o > a c . We now assume that each such fiber has a finite probability P(a,T) of failure at 
any stress a induced by a non-zero noise T: 



-i(^-l)], 0<a<a c 



P(a,T) = \ - 6XP L fW VJ' ----- i (5) 

{ 1, O > G c J 

As one can see, each fiber now has got a non- vanishing probability P(a, T) to fail under a 
stress a < a c at any non-zero noise parameter T. It may be noted that [unlike T in (2) or 



(4)] T is a dimensionless noise parameter. P(a,T) increases as T increases and for a > c c , 
P(cr, T) = 1. Unlike at T = 0, the bundle therefore fails at a < a c after a finite time r. Here 
we assume each fiber to have a fixed threshold a c , while their breaking probability at any a 
(< a c ) is due to noise- activated hopping over the barrier height (a c — cr). This differs from 
the earlier model studies [HI E] where the load distribution is noise induced. 

(a) Failure time 

At T 7^ and under any stress a (< a c ), some fibers fail due to noise and the load gets 
shared among the surviving fibers, which in turn enhances their stress value, inducing further 
failure. Denoting the fraction of fibers that remain intact at time t by U t , a discrete time 
recursion relation (see can be written as 



U 



t+i 



(6) 



where a/U t = L/(NU t ) is the redistributed load per fiber among the NU t surviving fibers 
at time t. In the continuum limit, we can write the above recursion relation in a differential 
form 



dU 
~dt 



a 

— exp 



T \o 



(7) 
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dU 



(9) 



for cr < a c . For a > a c , starting from U t = 1 at t = 0, one gets U t+ \ = from (5) and (6), 
giving r = 0. 

For small T and as cr — ► a Cl r ~ Texp (cr c /cr — 1)/T. This failure time r therefore 
approaches infinity as T — > 0. For a < cr c , one gets finite failure time r which decreases 
exponentially as a approaches a c or as T increases and r = for a > <j c . This last feature 
is absent in the earlier formulations [H]. However, all these features are very desirable and 
are in qualitative agreement with the recent experimental observations [Zj. This is also 
comparable with the phenomenological results from Griffith theory discussed in the earlier 
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section, although the crack size effect in the Griffith theory differs from that in the fiber 
bundle case. Our numerical study confirms the above analytic results [obtained using the 
continuum version of the recursion relation (6)] (see Fig. 1) well. 




Fig. 1. The simulation results showing variation of average failure time r against (a) stress a and (b) against 
noise T, for a bundle containing TV = 10 5 fibers. The theoretical results are shown by dotted and dashed 
lines [from eqn. (9)]. The insets show the simulation results for the variation of the fraction U of unbroken 
fibers with time t for different T values [1.2 (cross) and 1.0 (plus)] in (a) and a values [0.15 (cross) and 0.12 
(plus)] in (b). The dotted and dashed lines represent the theoretical results [eqns. (10 & 9)]. 



(b) Avalanche size distribution 

From the recursion relations (6) or (7), one can see that in each unit time interval a number 
of fibers break giving an avalanche size for the breaking. The avalanche size therefore is 
given by dU/dt and during the entire failure period r, different sizes of avalanches take 
place. Solving for U(t) from (7) one gets 



U(t) = — In 



T-t 



+ 1 



(10) 



LTexp(-l/T) 

employing the expression (9) for r. One can easily check that U(t) — 1 at t — and U(t) — 
at t — t (see Fig. 1). Also as t — > t c = r, U(t) decays as ln(r — t) ~ (r — t) 13 with (3 = 0+ 
from (10). Expressing dU/dt as the avalanche size m, one gets from (10) 

r-t 



m 



Texp(-1/T) 



+ l~r-t, 



(11) 



for T -> 0. 
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Here the avalanche size m can also be interpreted as the rate of breaking (dU / dt) and it 
varies with time as (r — 1)~ 7 , 7 = 1 as t — > £ c = r. Since r — £ corresponds to the cumulative 
probability D(m)dm of avalanches beyond t, one gets 

D(m) ~ m - a ;a = 2 (12) 

for the (differential) avalanche size distribution D{m). Also, the exponent of power law 
decay in eqn. (12) is independent of stress a and the noise level T which has been confirmed 
through numerical simulations (see Fig. 2). 




10" 3 10" 2 10~ 1 3 x10" 1 

m 



Fig. 2. The simulation results for the distribution D(m) of avalanches in the bundle with N — 1Q (averaged 

over 10 3 realisations): a = 0.2, f = 0.8 (triangle), a = 0.15, f = 0.8 (circle) and cr = 0.15, f = 1.0 
(square). The dashed line corresponds to a decay power 2.0. 

It may be mentioned that such avalanches manifest in the ultrasonic emissions during the 
propagation of fracture in the solid and the ultrasonic amplitudes are also observed to have 
similar power law distribution [T]. 
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IV. Simulation results for fatigue-failure in random fiber bundles 

In order to investigate the fatigue behavior in random fiber bundles we consider three different 
kinds of fiber strength distributions p(cr c ): (A) Uniform distribution of fiber strength where 
p(cr c ) = 1 for < cr c < 1 and p(cr c ) — for a c > 1, (B) Linearly increasing distribution 
of fiber strength where p(cr c ) = 2a c for < cr c < 1 and p(cr c ) = for o c > 1 and (C) 
Linearly decreasing distribution of fiber strength where p(cr c ) = 2(1 — a c ) for < a c < 1 
and p{(T c ) — for a c > 1. It has been already shown analytically [3|, from the dynamics of 
failure in all these three kinds of fiber bundles in the absence of any noise (vanishing T or T 



in (5)), the bundle's strength dc = 1/4 for model A, a~ c = y4/27 for model B and dc = 4/27 
for model C. We now consider the effect of the noise T inducing the failure probability 
P(a,T) = exp — i — lj for < a < cr c and 1 for a > a c , in the (fatigue) dynamics of 
such bundles, where a c is the strength of the individual fibers in the bundle. 

We have studied these numerically, using Monte Carlo method (for bundles having = 
10 5 or more fibers). We have considered bundles having the above three kinds (A, B and C) 
of p(cr c ) one by one. The noise induced failure mentioned above is realised only in a Monte 
Carlo way. Taking averages typically over 10 3 Monte Carlo runs the fraction of unbroken 
fibers U{t) at any time t at a fixed stress level cr(< 5^) is noted. At any a, the average failure 
time t (when U{t) = 0) is extracted. The form of the distributions and the variations of 
average time with noise T and stress a are shown for the three types of bundles. We find 
that r fits a form 

r = f«p(-i) exp(j| + i)-l (13) 

for all types of bundles (indicated by dotted lines in Fig. 3) . We find that this phenomeno- 
logical form (13) is indeed very close to the analytic result (9) for the fixed strength fiber 
bundle; it is somewhat approximate for these bundles and fits better for lower noise (T) and 
stress (a) levels. The avalanche size distributions in all these three models (A, B and C) 
have been studied numerically (see Fig. 4) and we find them to follow the same power law 
decay (12) with a ~ 2.0. 
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Fig. 3. Typical fiber strength distributions p(a c ) considered and the simulation results for fatigue behavior: 
(a) average failure time r vs. noise T (for three different stress values a) and (b) t vs. a (for three different 
noise values T) are shown for N = 10 5 fibers. The time variation of fraction of surviving fibers are shown 
in the insets for the three models: (A) with uniform p(a c ), (B) with linearly increasing p(<J c ) and (C) with 
linearly decreasing p(a c ); all having a cut off at a c = 1. The dotted lines in (a) and (b) corresponds to the fit 
with expression (13) where <r c ~ 0.245 in (A) (exact value=l/4 [3]), <j c ~ 0.370 in (B) (exact value=-\/4/27 
[3]), a c ~ 0.148 in (A) (exact value=4/27 [3]). 
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Fig. 4. The simulation results for the distributions D(m) of avalanches (m) in the three random fiber bundles 

with TV = 10 5 (averaged over 4 x 10 3 realisations): for model (A) with a = 0.07, T = 0.5 (square), for 

model (B) with a = 0.12, f = 0.4 (circle) and for model (C) with a = 0.04, f = 0.5 (triangle). The 
dashed line corresponds to a decay power 2.0. 

V. Summary and discussions 

First, we have studied analytically the macroscopic failure of a homogeneous fiber bundle 
model where each fiber has an unique strength (<r c ). At zero noise (T = 0) all the fibers of 
the bundle fail simultaneously for a > cr c , while at T ^ each fiber has got a non- vanishing 
failure probability [given by eqn. (5)] due to the thermal-like activation. The dynamics of 
failure of the bundle has been solved using the continuum version of the recursion relation (6) 
for global load sharing case. The resulting expression (8) for the average failure time (r) has 
qualitative features similar to that (4) obtained from the phenomenological nucleation rate 
theory applied for a Griffith's crack. Both the forms have got the desirable features that r 
decreases exponentially as a approaches a c from below and r ~ for a > a c . As mentioned 
already, although the above features agree qualitatively with the experimental observations, 
the precise mathematical forms we obtained here differ from the experimentally indicated 
forms [Zj. As time t approaches r, the fraction of unbroken fibers decay as (r — t) 13 , (3 = + 
and its rate of breaking grows as (r — t)~ 7 , with 7 = 1. The avalanche size distribution 
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D(m) is also obtained analytically for the dynamics. It is seen to have a robust power law 
governed decay behavior D(m) ~ mr a with a = 2. Our numerical results also confirm this 
behavior. Next, we have studied numerically the dynamics and the average breaking time 
t for bundles where the breaking strengths are not fixed and are given by the three simple 
distributions p(c c )- We find that for all the three cases, the average r fits well a form (13), 
which is very close to the analytic form for r in (9) for fixed failure threshold of the fibers. 
We have also investigated the avalanche size distributions in these models and obtained the 
same power law behavior, as for the fixed strength fibers. 

As mentioned already, here the noise parameter (T in (5)) can not be identified with 
temperature (T in (2)) which scales with the (crack) energy. In fact, although this failure 
model and its dynamics are applied here to classical breakdown phenomena occurring in the 
fiber bundle model or (classical) percolating solids fIJ, they seem to be applicable to quantum 
breakdown due to tunneling as well. Failures in quantum percolating solids beyond their 
linear conducting or insulating regime, has not been studied much (see however [8]). In fact, 
like the fuse (or dielectric breakdown) problems of percolating (or non-percolating) systems 
of conductor-insulator networks, one can think of the field induced breakdown of a quantum 
percolating system where the phase of the system is determined through two energy scales: 
Fermi energy e/ and the mobility edge e c . For e/ > e c the system is in conducting phase and 
it goes to insulating phase for e/ < e c . This metal-insulator transition at e/ = e c (in higher 
than two dimensional systems) and the scaling property of conductivity for e/ > e c have been 
studied extensively [HUH]. For the insulating phase (e/ < e c ), one can have electric field 
induced (Zener type) breakdown (similar to dielectric breakdown of non-percolating classical 
networks). This Zener breakdown of Anderson insulators or the quantum tunneling induced 
breakdown of impure (localised) insulators have not been studied much (see however |E1|12J). 
Unlike a (classical) fiber bundle model considered here, where all the fibers are in parallel, 
one can consider a dielectric composed of several elements in series having non-zero failure 
probability for each element due to quantum tunneling (like the noise-induced activation 
considered here). Any microscopic failure of such an element would result in increased 
field on the surviving elements and this in turn would enhance their failure probability. A 
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similar dielectric failure time (r) in such quantum or Anderson insulators is thus expected 
under electric field. Here a and a c would be replaced by e/ and e c respectively and T would 
correspond to the inverse tunneling length determined by the electric field (with the Planck's 
constant as the proportionality factor, incorporating the intrinsic noise) [HJ. 

Our study here for fatigue breakdown in the model fiber bundles shows that the aver- 
age failure time for the bundle at a stress value a less than the bundle strength a c (for 
homogeneous fiber bundle) or cT c (for random fiber bundles), above which the bundle fails 
immediately, decreases exponentially as a approaches o c or from below. This has already 
been observed in several experiments qualitatively. We have demonstrated this fatigue be- 
havior here both analytically and numerically for a homogeneous fiber bundle (section III) 
and also numerically for random fiber bundles with nontrivial strength distributions (section 
IV). We also believe that these observations will be useful in quantum breakdown phenom- 
ena. 
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